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Abstract 

The computation of PBH (primordial black hole) production from 
primordial perturbations has recently been improved by considering 
a more accurate relation between the primordial power spectrum and 
the PBH mass variance. We present here exact expressions which are 
valid for primordial spectra of arbitrary shape and which allow accu- 
rate numerical calculations. We then consider the possibility to have 
a significant part of dark matter in the form of PBHs produced by a 
primordial spectrum of inflationary origin possessing a characteristic 
scale. We show that in this model the relevant PBH mass is con- 
strained to lie in the range 5 x 10^^ S ^ M < 10^^ g. This is much less 
than the mass range coming from the QCD phase transition, allowing 
the two mechanisms to be easily distinguished. 



PACS Numbers: 04.62. +v, 98.80.Cq 



Introduction: A consistent paradigm seems to emerge in cosmology, of which dark 
matter is an essential ingredient. However, the nature of dark matter remains to 
be one of the most important open problems. It is expected to account for about 
one quarter of the present critical energy density, most of the remaining three 
quarters being dark energy whose nature is also unknown. There exist many 
candidates which are classified into cold, hot, and warm dark matter. Some of 
these candidates would signal new physics beyond the standard model of particle 
physics, such as relic neutralinos, the supersymmetric cold dark matter (CDM) 
candidate. Another possible CDM candidate could be primordial black holes 
(PBHs). The advantage with them is that their existence rests on known physics, 
general relativity and the presence of primordial fluctuations, independent of 
the mechanism which generates them. Indeed, it is well known that PBHs can 
be produced in the early universe due to the collapse of overdense regions on 
the linear size of the Hubble radius at the time when this corresponds to their 
Schwarzschild radius |I| . These fluctuations in the energy density must exist one 
way or the other in order to explain the origin of all inhomogeneities we see in the 
universe. Their existence, and in particular their primordial origin, is reflected in 
the presence of acoustic peaks in the cosmological microwave background (CMB) 
anisotropy, and is now well established by observations. Therefore, the generation 
of PBHs is unavoidable. Inflationary fluctuations are produced on a huge range 
of scales and the observation of PBHs formed after inflation could probe the 
fluctuations on scales where the CMB gives no information. It would therefore 
be complementary to the results which are extracted from the CMB data. 

PBHs were produced in the very early universe. PBHs generated before 
10~^^ s, corresponding to masses M < 5 x 10^^ g, have already evaporated 

by the present day due to Hawking radiation. PBHs with masses bigger than 
M* could, however, constitute a significant fraction, or even all, of the CDM. 
Unfortunately this possibility, though attractive, cannot be implemented with 
scale-free perturbations because they would lead to a negligible rate of PBH for- 
mation. However, in view of the naturalness of PBH generation, it is important 
to investigate whether other kinds of spectra could lead to a significant formation 
rate. 

General formalism: We will consider Gaussian primordial fluctuations, the kind of 
fluctuations expected in most inflationary models. The density contrast averaged 
over a sphere of radius R then reads 
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where 




computed using a top-hat window 



function. 
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P{k) is the primordial power spectrum, and WTuikR) is the Fourier transform 
of the top-hat window function divided by the probed volume Vw = |vri?^, 



WTH{kR) = ——^{sinkR - kR cos kR) . (3) 

The probability /3{Mh) that a region of comoving size R = ^^^\t^. has an averaged 
density contrast at horizon crossing tk in the range 6min < < ^max, is given by 



I3{Mh) = ^ / e d6 ^ Jl} e ^ (4) 

where cr'jjitk) = cr^(i?)|^^. Recently, the connection between the mass variance 
on very small scales and the normalization on present COBE scales has been 
improved yielding a substantially lower mass variance [@]. This formalism also 
allows for an accurate calculation of the mass variance for spectra with features. 
The quantity 5'jj{tk) = ^ P{k,tk) = k^ $^(A;,tfc) is often used; it is related 
to crjj{tk) for an arbitrary primordial power spectrum as ^ 

^Utk) = SUtk) T\kx, W^M dx , (5) 

where k^^'^{k,tk) = | F{k), tk -C teq (radiation-dominated phase). The time- 
independent quantity F{k) defines the primordial spectrum on super-Hubble ra- 
dius scales. The upper integration limit can safely be taken to be infinity due to 
the rapid convergence of the integral. The evolution of the fluctuations on small 
scales differs from the long wave evolution, and this difference is encoded in the 
transfer function T{k',t). At the present day we have = to, k = k^ = a^HQ 
and T{k',tk) 7^ T{k',to). This is why the relation between (Jnitk) and Sjjitk) is 
not the same as today. An accurate evaluation of the transfer function will be 
presented elsewhere IQ and gives 

T\kx,tk) = W^Hicsx) , = 

A pure scale-invariant spectrum corresponds to F{k) = constant, and a numerical 
estimate for a critical-density universe (flm,o = 1) gives @] 

alitk) = 5.37 4(4) = 6.63 5^(to) ■ (7) 

The prefactors in the second and last equality differ by a factor which takes 
into account the transition from radiation to matter domination. In particular, 
the constraint on a pure blue spectrum (spectral index n > 1) is weakened. The 
inclusion of a cosmological constant A with = 0.7 reduces the mass variance 
even further by about 15% Q]. This constant decrease is degenerate with an 
increase in 6min H) so we choose to make our numerical simulations for the 
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fiducial case f2m,,o = 1, ^min = 0.7. However, there is another important aspect 
of (^. For a spectrum with a feature at some scale ks, the shape of the power 
spectrum will translate in a non-trivial way into crfjltk) (see, e.g.. Fig. 1 below). 
Though the refinements described here make the prospects for the cosmological 
relevance of PBHs less promising since it reduces the probability to form PBHs, 
the situation looks bad only for scale-free spectra. Indeed, potentials with a 
feature on small scales could dramatically improve the PBH formation rate. This 
possibility was already considered in the past Q and we would like to take it 
up again in the light of the improved calculations done in ||^. 

PBH as a CDM candidate: It is not possible to have a significant fraction of the 
dark matter in the form of PBHs with a scale-free powerlaw primordial spectrum. 
We assume that the probability f3{M) represents the fraction of the energy density 
at the PBH formation time that is going to form PBHs with mass M. Their 
present contribution to the critical density, ^^pbh,o, is then given by 

npBHfiiM) h'' = Q.35 X 10'^ l3iM) {-jfj ■ (8) 

For a pure Harrison-Zel'dovich spectrum normalized to COBE, the mass vari- 
ance at Hubble-radius crossing, and therefore /3(M), is (nearly) constant on all 
mass scales. When we go to lower masses, the corresponding ^Ipbh,o increases 
cx M~5 because PBHs of lower mass form earlier in the radiation era, while their 
energy density decreases less rapidly than that of radiation. On the other hand, 
the COBE normalization of the spectrum results in (3{M) <^ 10"^'' so that this 
increase is by far insufficient and will result in a negligible contribution of PBHs 
to the present energy density. One way to circumvent the problem is to invoke 
a blue spectrum, n > 1, in which case n ^ 1.3 is needed in order to have a sig- 
nificant contribution with PBHs of mass M^, 0, §]. In such models, the spectral 
index is assumed to be constant over a huge range of scales. A more serious prob- 
lem is that as the Hubble mass at the end of infiation is lower than M*, a sharp 
constraint exists on the density of evaporated PBHs from their contribution to 
the 7-ray background, ^pbh,o ~ 10"^. Therefore, a blue spectrum that satisfies 
the bound from the 7-ray background necessarily yields flpBH,o ^ 1- The only 
way out is to boost the production of PBHs in a localized way. 

Microlensing gives no detection limits on the mass of massive compact halo 
objects (MACHOs) below 10^^ Mq. On the other hand, it is advantageous to 
form PBHs earlier, with masses not too far from M*. For larger masses, larger 
values for (3 are required in order to achieve the same flpBH,Q- We are thus looking 
for a spectrum which increases l3{M) in this region, but not in the mass range 
MpBH < M^,. This cannot be achieved neither with a blue spectrum, nor using 
a pure step-like spectrum with more power on small scales 0. In the following 
we shall present a model with a feature at some characteristic scale, where this 
is possible. 
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The model: We consider here a model based on an inflaton potential with a 
step in its first derivative. The adiabatic perturbations spectrum was derived 
analytically by Starobinsky ||^ and has been applied either to explain a possible 
bump in the matter power spectrum P{k) on a scale k ~ 125/i~^Mpc or on 



small scales as a cure to the dwarf galaxy problem |T^. In our case the feature will 
be placed on even much smaller scales, typically the scale corresponding to PBHs 
with mass M*. We expect that a deeper understanding of fundamental theories 
could give a distinguished scale from which the PBH formation rate could be 
calculated. In this model, the spectrum is defined by a characteristic scale kg (or 
mass Ms) and a parameter p, where is the ratio of the power spectrum on large 
scales with respect to that on small scales. In view of its rich structure in the 
vicinity of Mg, in particular its oscillations, one must use in this crucial part 
of the spectrum in order to compute aniM) correctly. The probability (3{Mh) 
to produce a PBH with mass M is obtained from the primordial spectrum in a 
non-trivial way according to 

a|(M) = F{bx) W^^icsx) W^^) dx , (9) 

with M = Muitk), MuitkJ = Mg, and tk is the Hubble radius ("horizon") cross- 
ing time k = {aH)t^. The non-trivial way in which the primordial spectrum F[k) 
translates into the mass variance aniM) is obvious from Fig. 1. 

We must distinguish three mass ranges: M <C Mp^ak, M ~ Mp^ak and M 3> 
Mpeak- For large masses M ^ Mp^ak, is unsignificant; in this range we have 

a (nearly) flat plateau normalized to COBE. For small masses, M <C Mp^ak, we 
are on the higher plateau satisfying ([TT|) (see below). Note that on both plateaus, 
corresponding to an effective spectral index n = 1, PBHs of arbitrary mass are 
equally probable in view of (H). The probability /9(M) that is obtained using 
(H,^ displays a pronounced bump around the mass scale Mpeak ~ 0.6 Mg, see 
Fig. 2. 

The crucial property is the drop on smaller scales M -C Mp^ak which, as we 
will show, is sufficient in order to satisfy the severe upper bounds set on the 
contribution of evaporated PBHs to the diffuse 7-ray background. 

Using the gravitational constraint (^), it is possible, for each value of the 
parameter Mg to find the value of p > pmin which yields flpBH,o < 0.3. This is 
shown in Fig. 3. Let us consider now the range of masses which have already 
evaporated by now and that contribute to the 7-ray background. Actually only a 
relatively small range AM with masses 2 x 10^3 g < M < 5 X 10^'' g contributes 



significantly to the observed diffuse 7-ray background ITT]. The fraction of energy 
density going into PBHs in this mass range obeys 

i t'PBHfl { i-Vlpeak ) \ ^v^peak / 
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Figure 1: The quantities crjiitk) (dotted line) and k^^'^{k,tk) oc 5'^{t}.) (solid line) 
are displayed for the same parameter value p = 0.5. The overall normalization 
is arbitrary (and different) for each curve. It is obvious that the two quantities 
crff{tk) and Sj^^tk) have very different shapes in the vicinity of the characteristic 
scale kg where cr'jjitk) exhibits a spiky structure. Their respective maxima have 
different locations. Note the same asymptotic ratio between large and small 
values of k (with tk <S teq). 



Taking flpBH,o (Mpeak) = 0.3, this leads to the constraint 



/ M \ ^ 

p. < 3.3 X 10-« — PiMpeak) . (11) 

V ^^J-peak J 

Note that we have used here the stronger constraint VLpBHfi{.M <C Mp^ak) < 10^® 
for each scale M <C Mpeak rather than the (integrated) constraint Qpbh o{^M) < 
10-8. 

A shift in the location of the characteristic mass Mg (and corresponding scale 
ks) to higher masses (still much smaller than 10"^ Mq) is allowed, but smaller 
values of p will then be needed in order to achieve flpBH,o = 0.3. However, this 
in turn increases the ratio -J^ which therefore gives an upper bound to Mpeak 

Ppeak 

around 10^^ g as can be seen on Fig. 3. Hence, requiring that ^pbh,o = 0.3, the 
7-ray background constraint can be satisfied for a range of characteristic masses 

M, (with M, ^ 1.6 Mpeak) 

lOM, <Ms< 10'' g . (12) 

For higher Mpeak, the 7-ray background constraint is no longer satisfied. An 
observational constraint on the allowed range of MACHO masses would therefore 
give a strong constraint on our model as it must correspond to the window (0). 

The bump for M ~ Mpeak is well localized, hence the PBH masses in this 
model indicate essentially the location of the characteristic scale Mg in the pri- 
mordial spectrum. Besides the various improvements we have used for the calcu- 
lation of PBH production, we note that the model presented here contains one 
singularity in the inflaton potential, in contrast to the more sophisticated version 
with two bumps considered in 0. 

It is interesting to compare these results with the approach in which PBHs of 
various masses M are produced by near critical collapse [|T2] at one single time, 
the horizon crossing time t^. corresponding to Mfj{tk) = Mpeak, with 

M = K{6- 6c)"' . (13) 



One finds for 7 = 0.35, following []T3|, 



dM Mpeak \ Mpeak J ^ \ Mpeak J ^ ^ ' 

This is turn yields 

^PBH^tk^^J = 0.8 l3{Mpeak) , (15) 

a result analogous to that obtained when PBHs of various mass M are pro- 
duced at different times tk with M = Mnitk)', in the latter case one obtains 
^PBH,tot{tkp,,k) = P{Mpeak)- Hcucc, our Constraint on ^PBH,o{Mpeak) using (D 
is stronger than the requirement ^pBH,tot{h) < 0.3 based on (pISl). Note that 
Mpeak-Mn in the notation of is the Hubble mass corresponding to the peak 
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in the quantity /?(M), not in the primordial spectrum itself, due to the distinc- 
tion one has to make between the quantities anitk) and Snitk)- We stress also 
that the above mentioned approach does not take into account the properties 
of the feature in the primordial spectrum, in particular the width of the corre- 
sponding bump in /3, besides the value (3{Mpf,ak) which serves as a kind of overall 
normalization of the PBH abundance. On the other hand, it does account in a 
consistent way for the production of PBHs also in the monotonically increasing 
part of /3{M) around Mp^ak, M < Mp^ak- Obviously, the more spiky the pri- 
mordial spectrum, the better this approach is expected to be, and our model 
provides such a concrete spectrum based on an underlying inflationary dynamics. 
An accurate determination of the mass spectrum would therefore yield valuable 
information, not only on the underlying primordial perturbations spectrum but 
also on PBH formation itself. 

In conclusion, we have considered a power spectrum of inflationary origin in 
which PBHs constitute a significant part of CDM. In fact, it has already been 
speculated that the QCD phase transition could lead to PBHs of about one solar 
mass [Q, in accordance with microlensing observations made some years ago. In 
our model PBHs are generated with much smaller mass. Our analysis is based on 
an accurate calculation of the mass variance at the PBH formation time tk and 
of the corresponding (3{M). Though it was found recently that a more accurate 
computation of the probability to produce PBHs actually leads to a decrease of 
this probability, this letter shows that one can still have a significant Qpbh,o by 
using primordial spectra in concrete models with a characteristic scale on very 
small scales, scales much smaller than those usually probed by CMB anisotropy 
or large scale structure formation. 
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Figure 2: The quantity /3(M) is shown for our model containing a jump in the 
inflaton potential derivative for the parameters p — 7.913 x 10""^, Mg = lO^^g. As 
can be seen, (3{M) acquires a well localized bump in the vicinity of Mg = M{tk,). 
Note that wc have Ms 1.6 Mp^^k and = (3{M < Mp^ak) ~ 10"^^ x (3{Mp^ak), 
a value sufficient to avoid the severe constraint from the contribution to the 7-ray 
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Figure 3: The allowed region in parameter space {p, Mg) is shown. The solid hne 
indicates the points for which ^pBHfl{Mpeak) = 0.3. Below the solid line, the 
gravitational constraint at Mp^ak is violated and this region is therefore excluded. 
Below the dotted line, the 7-ray background constraint is violated. It is seen 
that for lO^^g < M^, the allowed parameter values p yield ^pBH,o{Mpeak) < 0.3 
which becomes rapidly negligible with growing M5. For the values of p shown 
here, Ms 1.6 Mpeak- 
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